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Paolo AmorcQ 

Facultad de Ciencias, Universidad de Colima, 

Bernal Diaz del Castillo 340, Colima, Colima, 

Mexico. 

A method devised by the author is used to calculate analytical expressions for one-loop integrals 
at finite temperature. A non-perturbative regularization of the integrals is performed, yielding 
expressions of non-polynomial nature. A comparison with previously published results is presented 
and the advantages of the present technique are discussed. 

PACS numbers: 


I. INTRODUCTION 


In this paper we present a new technique which allows 
to evaluate non-perturbatively one-loop integrals occur¬ 
ring in finite temperature problems in quantum field the¬ 
ory. Such integrals typically involve series, due to sum¬ 
mation over the the Matsubara frequencies, = 27rnT, 
with n = 0,±l,... P|, which cannot be done analyti¬ 
cally. 

The method that we propose here is derived from the 
powerful ideas of the Linear Delta Expansion (LD^ and 
of Variational Perturbation Theory (VPT) (see [J and 
references therein); we briefly sketch how such methods 
work. 

First one interpolates a non-perturbative lagrangian 
(hamiltonian) £ with a solvable one £o, which depends 
on one or more arbitrary parameters; the interpolated la- 
gragian is then split into a “leading” term, chosen to be 
Co itself and into a “perturbation ”, C — Cq. Although the 
“perturbative” term will not be a-priori small, a pertur¬ 
bative expansion is carried out and physical quantities 
are calculated. To finite order in perturbation theory 
one thus obtains expression which contain an artificial 
dependence upon the arbitrary parameters: however if 
the expansion were to be carried out to all orders such 
dependence would have to cancel out provided that the 
perturbative series converges. In order to minimize such 
dependence one applies the Principle of Minimal Sensi¬ 
tivity (PMS) 0 , which selects the value of the arbitrary 
parameter for which the physical quantity is less sensi¬ 
tive to changes in the parameter itself. The parameter 
determined in this way is normally such to make C — Co 
a perturbation and depends on the “natural” parameters 
present in C. For this reason, the expansion obtained by 
using this method does not provide, to a finite order, a 
polynomial in the parameters of the model, as would be 
the case for a genuinely perturbative technique. 

The ideas behind the LDE and the VPT are really very 
general and have been successful in dealing with a vari¬ 
ety of problems of quite different nature, ranging from 


quantum field theory to classical and quantum mechan¬ 
ics. In this paper, in particular, we pursue the application 
of some results which were recently obtained by the au¬ 
thor 0: in that paper the author developed a method 
to accelerate the convergence of certain class of mathe¬ 
matical series (including the Riemann and Hurwitz zeta 
functions) and proved that the new series obtained in 
such a way converge exponentially to the correct results. 
As we will see in this paper, some of the series treated in 
0 are relevant for the evaluation of one-loop integrals at 
finite temperature and will be used to obtain arbitrarily 
precise analytical approximations to such integrals which 
are not polynomials in the inverse temperature. 

The paper is organized as follows: in sectionwe re¬ 
view some of the results of 0 and explain the method; 
in section llni we apply the results of section|n]to the cal¬ 
culation of one-loop integrals at finite temperature and 
compare them with the results in the literature [1,IE00; 
finally in section Hvl we draw our conclusions and discuss 
further possible applications of this method. 


II. THE METHOD 


In this section we describe the method of 0, which 
allows one to accelerate the convergence of certain series. 
One of the applications which was discussed in that paper 
is to the zeta-like function defined as 


n=0 


1 


( 1 ) 


where s, u and ^ are real parameters. The special cases 
obtained by taking {u,^) = (1,1) or u = 1 correspond 
to the Riemann and Hurwitz zeta functions respectively. 
Further examples of application of this method are con¬ 
tained in 0 and the interested reader should refer to that 
paper for further reading. 

We write 
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and by simple algebra we convert it to the form 
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having introduced the definition 


A(n) 


1 + A2 
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We stress that eq. © is still an identity and that no 
approximation has been made so far. A here is an ar¬ 
bitrary parameter which was introduced “ad hoc”: this 


procedure is typical of the Linear Delta Expansion (LDE) 
and of Variational Perturbation Technique (VPT) @]. 


Clearly, when the condition |A(n)| < 1 is met, one can 
expand eq. m in powers of A; such a condition requires 


that 
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< 1, i.e. that A^ > 


We can therefore write 
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and convert the series © to the equivalent series: 
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By interchanging the sums and performing the sum over n we finally obtain the expression 
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We stress that eq. 0 is still exact for A^ > and 
that it converges exponentially. Actually this equation 
describes an entire family of series, all converging to the 
same function, although with different rates of conver¬ 
gence. When the sum in eq. 0 is truncated to a finite 
value, a fictitious dependence on the arbitrary parame¬ 
ter A is generated: in order to minimize such dependence 
and to obtain a series with an optimal rate of conver¬ 
gence one can apply the Principle of Minimal Sensitivity 
(PMS) 0, corresponding to selecting the value of A for 
which the derivative of the partial sum with respect to A 
vanishes. 


III. APPLICATIONS: ONE-LOOP INTEGRALS 
AT FINITE TEMPERATURE 

We will follow the notation set in @ and write the 
expression for a general one-loop integral at finite tem¬ 
perature in the form 


J(m,a, 6) = ^ 


n = —oo 

n y 0 


d^k 

(2^ [fc2+a;2+m2]‘ 


To leading order one obtains the result 


being the scale brought in by dimensional regulariza¬ 
tion, a and b being integers (a > 0 and 6 > 0). Following 
we also define 
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which can then be used inside eq. 0 . The precision of 
these formulas has been investigated in Q , showing that 
extremely precise results are in general obtained already 
working to low orders. 


and obtain 
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where 


SiM,l) = Y, 


{n^ + M^y 


( 10 ) 


and I = b — a — D/2 and D = 3 — 2 e. Depending upon 
the value of I, the series of eq. mu could be divergent 
and therefore need regularization. 

Eq. lion is clearly of the form considered in eq. and 
one can write 

SiMJ) = +a2,1, M^) , (11) 

where 


({2,1, M^) = 


and 

^fe(A,2,Z,M2) = 


M2* 
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Once more we stress that although the series all con¬ 
verge to the same result independently of A (provided 
that A > ^^), the partial sums, obtained by truncat¬ 
ing the series to a finite order will necessarily display a 
dependence on the parameter. Such dependence, which 
is an artifact of working to a finite order, will also make 
the rate of convergence of the different elements of the 
family A-dependent. Since it is desirable to obtain the 
most precise results with the least effort, one will select 
the optimally convergent series by fixing A through the 
“principle of minimal sensitivity” (PMS) |^. The solu¬ 
tions obtained by applying this simple criterion display 
in general the highest convergence rate and, once plugged 
back in the original series, provide a non-polynomial ex¬ 
pression in the natural parameters. 

In the present case, the “natural” parameter in eq. m 
is M^ and to first-order the PMS yields 


'^PMS 


C(2(l + 0) 
C(20 


(13) 





FIG. 1; C(2, 3/2,^) calculated numerically (solid line), with 
the PMS to first order (dashed line) and using A = 0 (dotted 
line). 


In fig. n we have compared C(2, 3/2,^) calculated nu¬ 
merically with the first-order approximations obtained 
by using eq. m with A = X^p'lfg and A = 0, which are 

d'L(2.3A{) = F7r+(f(3)+ff,5))V2 (15) 
Cfio(2,3/2,e) = ^ + C(3)-^?C(5) , (16) 

where the second and third term in the second equation 
can be obtained by Taylor expanding the first equation. 

In table^we have calculated ({2, 3/2,1) using eq. itT^ 
to order N with X^p'j^g (second column) and with A = 0 
(third column). Our formula, using the optimal param¬ 
eter obtained to first order, converges exponentially to 
the exact value, whereas the formula corresponding to 
A = 0 is actually useless (that formula is indeed limited 
to ^ 1). Notice that ^(2,3/2,1) calculated with the 

first 10® terms in eq. CQl gives the result 1.5124349215 0: 
the same precision is reached with our improved series 
with only 27 terms. We believe that this result by itself 
is sufficient to illustrate the strength of our method. 

As an example of the application of our formula, we 
consider the case studied in eq. (12) of @], i.e. 


On the other hand, if the value A = 0 is chosen, then 
one obtains the series 


({2,1, M^) 
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which, to a finite order yields a polynomial in M. Notice 
that such series corresponds to the expansion used in . 
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which is essentially the one-loop self energy. 

The first integral is finite and equal to —MT/{Air) 
whereas the second term can be written as: 
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TABLE I: C(2,3/2,1) calculated using eq. 11211 to order N with (second column) and with A = 0 (third column). The 

underlined digits are correct. 

JV _Cpms(2,3/2,1)Ca=o(2, 3/2,1) 

1 1.4728636067646540152245808470808148333385501418201 0.64666527044453939590268993182589873867936491308763 

10 1.5124 141548516402130603869952574426525740627123155 3.4054827271508252417934796716746660253808017949304 

100 1.5124349215502030648030954892350003 557973564160761 7.189509053068942853725754726982871981048977697327 

200 1.5124349215502030648030954892350003618366746247481 9.5161848472956355375285584079161971312864869056589 


J{M, 0,1) = 2-2-2'= p^_j^2 + e) 

I- 


+ C(2, -1/2 + e, m2) 


(18) 


where the divergencies are contained in C(2,—1/2 + 

e,M2). 

By using our eq. and retaining only the divergent 
terms and those independent of e we obtain 


, mT T2 m2 

I{m) = -1- - 

47r 12 IGtt^ 


1 , 47rT2 

- + 7 - log 
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'C(3) 


‘C(5) 


5m®C(7) 

8(27r)4r2 “ 1024 7r6 32768 tt^ 

7mi°C(9) 


(19) 


262144 MOTS 

which reproduces the results of 0, which, however, were 

I 


considered only up to order m^. As anticipated, the for¬ 
mula obtained is a polynomial in m. 

We will now use the optimal series 03 to improve this 
result. We first notice that since 


dC{u,s,^) 


= -s C{u,s + 1,0 


( 20 ) 


a divergent series can be related to a convergent series by 
taking repeated derivatives with respect to the parameter 
0 Indeed by applying eq. twice we can write the 
general expression 


J{M,a,b) = 


2l+2“-2f'+2'= ^|+2a-2b-He y4+2a-2 6p(-3 - a + b + c) 


r(6)r(|-e) 


1 “h a — b -\- 


3-2e 


—a b -\- 


-3 + 2e 


I J dM"^ j dM'^ I^-M-2 (2-a+b+^^ ) ^(2^ l_a + b + e, M^)'^ | (21) 


where the last / function is now fully convergent when divergent contributions. 


a = 0 and 5=1. Since this formula requires a double 


By applying eq. 03 to first order and using the op- 


integration of a convergent series, the result is determined timal value of A given by eq. 03 we obtain the simple 
up to two constants of integration, which will contain the result 


Jpms{M,0,1) = T2 -Ki M2- 


C(3)^/^ VC(3) + M2C(5) ^ 

C(5)^ 


K2 


( 22 ) 


where ki .2 are the constants of integration independent 


of M. In order to determine these constants we Taylor 
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expand this expression in M and then use M = m/{2'KT) 
thus obtaining 


J(M, 0,1) 


47r2 128 7r4T2 (-(5)2 

m? (,{2,)^ 7ti®C( 5) 5 m® ((5)^ 

87r2C(5) ~ 1024 tt^ 32768 tt® C(3) 


7mi0C(5)^ 

262144 7rior®C(3)^ 


(23) 


Eq. ^ is a quite remarkable formula: indeed, al¬ 
though it has been derived by applying our method to 
first order, it reproduces correctly the terms going as m^ 
and m® and it also provides the coefficients of the higher 
order terms apart from factors depending on the Rie- 
mann C, function evaluated at odd integer values Of 
course, we can let the perturbative result guide us fur¬ 
ther and use the fact that lim ^(n) = 1 to write the 

n—♦■oo 

improved formula 


This expression can be now compared with the pertur¬ 
bative result and the constants of integration can thus be 
extracted 


Ki 


K2 


- + 7 - log 
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1 , C(3)^ 
12 C(5)^ 


47rT2\ 
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2C(3)^ 

C(5) 
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^PMS2(A7, 0,1) = 


13^2 

12 4 


T2 

l6 


- + 7 - log 
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47rT" 


- 2 


(l6 - 2 (-1 -F C(3)) + M^ {-1 + 


(24) 


where the notation PMS2 has been introduced to dis- By putting the pieces together we finally obtain our 
tiguish it from the previous formula. Notice that both approximation to /(m), given by 
equations (El and dMI) are non-polynomial in m. 
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IpMS2{M) 


MT^ 
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m 2 r2 
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- + 7 - log 
e 




T2 

16 




2 (-1 + C(3)) + M® (-1 + C(5))) 


(25) 


When negative values of m^ are considered, correspond¬ 
ing to a spontaneously broken phase, these results hold 
for T > To = —m2/(47r2). 

As a further application of the method described in 
this paper we consider the integral 

= [ ^ 2 TTY- 2 e log (</>)) (26) 


where fco = 27rnT, for all n integers. The effective poten¬ 
tial to one-loop is expressed in terms of this integral. 
Given the relation 


I{m) 


1 

m dm 


(27) 


it is not necessary to calculate J{m) which can be ob¬ 
tained as 


J{m) 


= y I(m) m dm + p , 


(28) 


where p is constant of integration independent of m. 

By using IpmS 2 of eq. El we obtain the simple ex- 
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pression: 


JpMS2{'m) = p + 

T 

Gtt 




2 rT^2 

m 1 — 


m^T 

12 TT 


(m2 + 4 7r2T2)"/' 


6471^ 


- + 7 - log 
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47rT" 




- 2 


768 tt^T^ 
^8 


(C(3)-l) 


m 


(C(5)-l) 


8192 TT^r^ 

which can be expanded in powers of m to give 

dTT^r^ TO2y2 


(29) 
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6477^ 
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? 24 

i + 7 - log 
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12 TT 
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®C(3) 


^C(5) 


768 7r4 T2 

j^io 

65536 TT® re 


8192 TfS r4 


(30) 


which reproduces exactly eq. (3.7) of 0, provided that 
p = ^ —h constant. This last constant is independent 
both of m and of T. 

We are now in a position to compare our result of 
eq. 12911 with the high temperature expansion of eq.(42) 
of ref. 0 which reads^ 
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2VT^ 


TT 

1677^ 
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127r 
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y 4 

1) C(2fc + i) 


fe=i 


r(A: + 1 ) r(fc- 
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(31) 


after setting the chemical potential to zero^. We have 
defined y = m/T. 

As we can see from table HD the result of 0 for I{m) 
would essentially provide the expansion of eq. JSU and 
thus reproduce the results obtained long time ago by 
Haber and Weldon 0 . We regard this procedure as “per¬ 
turbative” meaning that, when the series in eq. eu is 
truncated a polynomial in powers of m/T is obtained. 
On the other hand, our simple formula, eq. 129|l repro¬ 
duces correctly the perturbative expansion of [g up to 
order (m/T)®, and up to a factor involving the C, func¬ 
tion calculated at odd integer values, which however tend 
to 1 for large values. This makes our simple formula quite 


TABLE II: Comparison between eq. (EU and the simple ap¬ 
proximation of eq. (EU. 


N 


eq. 1311 


JpMS2/T''‘ 


T 

1 10 


- T J 

:f] 


C(3) 

768 7r« 
C(5) 

8192 -n-B 
C(7) 
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C(5) 
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3145728 
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FIG. 2: Comparison between JpmsiIT'^ and eq. 131 1 taking 
t^max — 10 (dashed), kmax — 20 (dotted) and kmax — 30 
(dot-dashed) terms. 


precise. We remind the reader that eq. (P|) was obtained 
by applying our method to Erst order and that, given the 
convergence of the method, drastic improvements in the 
quality of the approximation are expected if higher orders 
would to be taken into account. Although the calculation 
of higher orders with our method would be an interesting 
issue by itself, we will leave it for future work. 

In fig. El we have plotted JpmS 2 /T'^, defined by taking 
out the terms up to order {m/T)'^, and the similar func¬ 
tions obtained from eq. m, performing the sum to three 
different orders {kmax = 10, 20, 30). The horizontal scale 
isy = m/T. We leave to the reader any judgement on the 
quality of our approximation. Notice that the series of 
Haber and Weldon should correspond to our variational 
series for A = 0; as we have seen before the criterion of 


convergence for a series is that > 


M^-l 


which for 


A = 0 can be fullfilled only if m/T < 27r. Looking at the 
figure we indeed see that when we get close to this value 
the sum of 0 becomes ill-behaved. 


IV. CONCLUSIONS 


2 We divide by an overall factor of 2. 

® In such limit the hypergeometric functions in the original formula 


of 0 all go to 1. 


The method that we have described in this paper is 
quite general and probably could be extended in the fu¬ 
ture to deal with a larger class of problems than the one 
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presented here. We have proved that, by using varia¬ 
tional techniques it is possible to estimate analytically 
and to an arbitrary degree of precision series which are 
difficult to evaluate with standard techniques. In partic¬ 
ular, we have shown that by relating divergent series to 
convergent ones through repeated derivatives and then 
matching the divergences contained in the constants of 
integration with the ones coming out of the perturba¬ 
tive calculation, it is possible to construct a really non- 
perturbative regularization. Such regularized expressions 
are very accurate even at low orders and low tempera¬ 
tures. There is a huge literature dealing with field the¬ 
oretical problems at finite temperature and we feel that 
this paper can provide a quite general and useful tool 
to attack many of these problems. We also believe, 
although at present is still to be confirmed, that the 
method that we have described could be useful in the 
non-perturbative calculation of the Casimir effect, for 
which zeta function regularization is a well-established 


technique (see for example i). It would be quite in¬ 
teresting to see if non-perturbative expressions for the 
Casimir effect could be calculated analytically. We hope 
to apply these ideas to such a problem in the near fu¬ 
ture. As a final remark, we like to stress that despite 
the simplicity of the ideas that we have illustrated, all 
the results that we have obtained are fully analytical and 
improvable to the desired level of accuracy. 
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